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ABSTRACT 


This report is concerned with the approximation of acoustical noise 
fields by the use of modulated pure tone sources. Given a specific random 
signal the methods of orthogonal expansion and Fourier expansion are studied. 
The integral equation arising from the orthogonal expansion is solved for 
the normal functions when: a) the power spectrum is given as the ratio of 
two polynomials, b) Markoff Autocorrelation function, c) stationary band 
limited Gaussian white noise of mean zero, d) band limited stationary normal 
white noise on a finite interval, e) white noise. Deterministic and random 
frequency and amplitude modulation and random switching are examined as 
ways to broaden spectral components. 
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I. Introduction 


Because of the many applications, it is found useful 
to represent a random process In a suitable series of 
particular functions over seme finite or infinite interval. 
Considerable study has been devoted to the problem of 
expansions (Ref. 1,2,3)* Orthogonal expansions have 
received particular attention (Ref. 2). 


The present report deals with seme aspects of the 
expansion of random functions into orthogonal series. 
Consideration is also given to how to broaden the spectrum 
of a single-tone signal by modulation of various sorts. 
Some attention is also given to optimal approximations of 
various statistical characteristics of actual random noise 
fields. 


Manuscript released by author October 7* 1964, for publication as 
an RTD Technical Documentary Report. 
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II* Representation of Acoustic Noise Fields 


a) Orthogonal Expansion of Random Functions 

It Is desired to represent a randomly varying sound 
pressure p(t) (see Figure 1) as a series of orthogonal 
functions with orthogonal random variables as coefficients. 



I 



t 


Figure 1 , Typical Recording of Random Sound Field 


In particular, we require 




1 n ■ k 
0 n * k 



(t)dt 





To achieve these results, consider 

P(t) pV)* n (f) -I IX.I 1^1 ♦ < (t)t*(t')t n (t') X.X. 

8 ,k 


2 














Integrating both sides of this result from t- 0 to T and using 
the second of (2), we find 


jTVutV K iv )dt ' 



-I i\i 


X X 

0 n 


TV>Hng expect ion and lf ** n g the first of (2) gives 


J R(t,t')* (t')dt' - |X | 2 ♦ (t) 

n n n 


R(t,f) - K{p(t) p(t')J 


o< t < T 


(» 


Thus we have such an expansion as (2) provided ♦ (t) 
solutions of (3) and X is the eigenvalue associated with 

*.<*>• 

Fran (1) and (2) it follows that 


X n * l^ri V *> *<*)« <*> 

Therefore if we knew the eigenfunctions t^( t), and the 
eigenvalues X , and possessed a realised value of p(t), 
we could calculate the realised values of the orthogonal 
random variables X . All that is needed to find X and 
t (t) is the autocorrelation function R(t,t'). Bi?s is 
either given or is determined from an actual recording 
of p(t). 


If p(t) is a gauss ion random process, the ra ndom 
variables X will not only be orthogonal but will be 
Independent random variables with a normal probability 
distribution. The problem of determining the probability 
distribution function for p(t) when those of are given 
is an extremely difficult one and has been solved only in 
the case of gausslan variates. It is to be noted from (4) 
that if p(t) is an Independent process, X will be a 
normal variate. One can easily prove this by application 
of the central limit theorem of probability theory. 
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b) Fourier Expansion of Functions 

Because of the difficulty of solving the Integral equation 
(3) and because of simplicity, random functions are expanded 
into Fourier series in most practical situations. Unfortunately, 
the random coefficients are not orthogonal on (0,T) except for 
periodic processes. We can show, however, that the coefficients 

are asymptotically orthogonal as the length of the recording 
increases. 


Consider a stationary random function p(t) defined on the 

interval (0,T). The Fourier expansion of this process on (0,T) 

is 

m 

p(t) cos « k t + b k sin w^t), - ku (5) 

k-0 



O o 

T 

§ J p(t) cos kW Q t dt,k -1,2,... 

o 

T 

§/ p(t> sin ftt,k -1,2,... 



Now 


E **k a j* “ i? ip( t )p(*0 cos ka t cos iw t'dtdt'} 






cos kw t cos iw t'dtdt* 
o o 




fofo 


cos ku t cos Uti t'dtdt' 
o o 


T - t-t* 

Changing variables the double integral beccoes 

E{ Vi ) " &JI cos J*~ X R(t) COS (t>fT ) dt 

•t 

Let A - t/T, then 
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TE 


‘Vi 1 ’ 


i Td-w 

kl cos 2rricXdX / R(t) cos 2ri(~ * dT 


At 






T(l-X) 


■ cos 2vkXaX (cos 2 ttJX / R(t) cos 

p T(l-X) 

- sin 2rriX / R(t) sin dr) 


2 tt.It 


dr 


2 ttJt 


As T X ^ 0 


cos 2irkXdX (cos 2rrJX* j- S(m^) + 


+ sin 2irJX ; (0) } 


(7) 


s <“i> 


2wl 


■ Spectral density of p(t) and ti w . 
unless i Is chosen such that 11a w«-» w > 0. 

i,T -» «» 


(7) nay he evaluated to give 


TEte^} -♦ S(“,) 6 lk 


The noznalized a ■ 



obey 


0 81 T -> » 


( 8 ) 


E{ Vi , iTr 5 ik ' k ^ 1 


(9) 


In like wanner one can shoe that 


E 


<\V 


E(S S } -♦ 8 

1 kJ ' ki 


•‘Vi’-'oi 

E *oV- # 01 


k,i > 1 


( 10 ) 


Thus the expansion coefficients are asymptotically orthogonal. 


5 




It is of value to compare various Fourier expansions for 
different lengths of T of recordings of p(t) with the exact 
orthogonal expansion on (0,T). This would give some idea of 

the length of record to he taken so that a Fourier expansion 
would he as good as the exact one. 



Expansion of a Function. 


a) General Results for Processes with Rational Spectr* 1 
Density Functions 

We confine ourselves to stationary random sound fields whose 
spectral density functions are rational functions of frequency. 
In this case the Integral equation becomes 


-.T 

J R(t-t') t (f)df - |* | 2 t (t), 0 < t < T 



By the Wiener-Khintchine theorem, we know that 

R(T) m »f S(w) elTW<3w 

S(w) m J' r( t ) e iTW dr 


( 12 ) 


~>(w) is the power spectral density of the random sound pressure 

p(t). We now suppose that S(w) is a rational function of w so 

that it can he written as the ratio of two polynomials in w, 
namely. 


3(w) . UiMl 

D ((iw) 2 ) 



D is at least one degree greater than V. The argument of the 
N and D pol ynomi a ls was chosen as a square of w because of the 
general requirements placed on S(«) as a power spectral density 
function. We wish to show that when the process has such a 
spectral density as ( 13 ), the functions *r (t) satisfy a Unemr 
constant coefficient differential equation. 


We start hy inserting (13) into (3) to obtain 
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bf f e l(t ’ t,)h V - |X | 2 t (t) 

J o D((l») ) n n n 


(Ik) 


Let 


H((l«) 2 ) 


q/2 

I \ -(lB)2k 


K(i») 2 ) 


p/2 p > q 

2 k 


lb ^ ' liu) 


(15) 


a and b are constants. Bov 
with the operator 


operate upon both sides of (lk) 


P/2 


D( f?> -I 


kmO 


2 k 

to give 
dt 


4 D((ia) ) “ 


(16) 


- i\i*«!?> ♦„<*> 

2 p°° 

Multiplying by H((lw) ) under the Integral j is equivalent to 
operating on the integral vith N(^r). "*** Thus (l6) can be 

written as 


But 


f 0 [ "‘3? ) /'" i(t ’ t ) “ d “ 1 *n (t ' ,dt ' 

- I\l* »4V*) 

J e 1 ( t *'t _ 27T B(t-t') 
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where 6(t-t') is the Dirac delta function. Using (l8) in (17), 
ve find 

2 T 2 

»(J*(t-v) * n (f)dt- - |X n | 2 Dffp) t n (t) 


or 



is the basic differential equation that must be satisfied 
by the expansion functions ♦ (t). 


The solution of (19) will contain A and p arbitrary 

constants A . Khan this solution is substituted back in the 

integral equation (11), it will be seen that the integral equation 

cannot be satisfied except for certain values of A, the 

characteristic values, and the constants A must satisfy certain 

conditions. These conditions show that for each n there is only 

one independent constant. Thi« one constant is to be determined 

from the normalizing conditions on t (t), namely, 

& 


f * w ♦*<t)dt. ( ; i; j (20) 

o 


Thft procedure Just outlined will now be carried out explicitly. 
The differential equation (19) becomes 


q/2 

I 1 

k-0 


2 k 

4 V*> 


dt 


2 k 


P/2 

Ts 


2 k , , 
d f n (t) 


k-0 


dt 


2 k 


( 21 ) 


Try a solution 

if (t) - A e^n* (22) 

n n 

A is a constant and a is a parameter. Putting (22) in (21) 

we find that fl must b8 a root of 

n 
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(25) 


q /2 

I 

kvO 


2 k 

b a 

k n 


P/2 

|X /I vfi 

k»0 


2 k 


Ihis equation will have, in general (special situations 
taken care of), p roots £ Q - ♦ 0 (X ), k - l,..,p/ 2 . 

of the roots are simply the negatives of the other half, 
general solution of ( 21 ) is 


be 

Hair 

A 


t n (t) 


P/2 

l 

knl 


0 , t - -Q t 

(A Ok * + V ' “ > 


(24) 


The autocorrelation function is given by 


R(t) 




H((iw) ) ixw. 
— i ~T— e du 

D((i«) ) 


(85) 


This integral can be evaluated by contour integration using 
semicircular arcs and yields 


R(t) 


P/2 

I 

k -1 

P/2 


k 



r, T 

Be k 
k 


T - t-t'X) 


(26) 


T < 0 


w. are the roots of 


D((iw) ) - 0 


(27) 


lying in the upper half of the complex plane as shown in 
Figure 2 
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There are p roots of (27) + <*» k > 4 “ !>*••> p/ 2 * ® ie 

B. are given by 
k 


B 


N((iu» k ) ) 


k 2 

dD((iM) ) 


dw 


(28) 


U) - w. 


Inserting (24) and (26) into (11), ve obtain 


P/ 2 


P/ 2 



iu (t- 
B e s' 

'8 


f) V fA + 0 V tf 

k 4 lA nk * ^ 


+ A e ) dt' ♦ 

nk 



P/2 


-iw (t-t 
Be s 




.W aV.-^' 

ok nk 


) dt‘ 


P/2 


l\ 


I e Vt 


- -Q , t. 

+ A e nk ) 
nk 
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Carrying out the Integrations ve obtain 



+ e 0 *** 


B A 
i nk Q 


°nk V V “ e” 0 *^ 

( .-1 u * A, A B « A nk -fl -1 m 
nk s T ~k nk s 



B A 


+ e 


Qnk 


s nk ft + 1m 
nk s 



B A 


ft -fl n k t 


8 nk -Q +iw 
nk 8 


i £ v; 


lM t 

e 8 

nk 0 , - 1 m 
nk 8 




B A 


- e 


ifc^t 


8 nk -fl , -lu 
nk i 


v V ♦ e (°nk +i ‘^ T _ifc, t 

+) BA, ---- e s 

L A 8 nk 0 , + 1 m 
8 k nk s 


- e 


(“°nk + 


+ II B . A »k '“-v ♦ *■, 


-lM t 

e 8 


|x /£ K* #%lt + • 


This equation can be satisfied identically In t only If the 
sum of the coefficients of the different functions of tine 
are equated to zero. This gives the system 


(a) 


P/2 .+ 

y, ( .k. 

A s l o -iw 

*51 nk 8 


nk 


-0 ■-1m 

nk s 


} ■ 0| 8 ■ 1>•••>p/2 


0 >) 


P/ 2 A + e (0nk ♦ ^s) 1 

V _ , nk 



B (• 

8 


a, + in 
nk s 


A e 
nk 


(-link ♦ 


-0 . +1» 
nk s 


) - 0 


P/2 


(c) 


I V 

a«l nk 


8 ■ 1|...jp/2 

i«- - JrW 1 B . ’ M 2 V 

s nk s 


k ■ 1 ,...,p /2 


(29) 
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1 


1 


P/2 


I 


A*. {■ 


_ - ** 0 , - 1 **. 

»mO ok S 


-fi 


nk 


7ir ,B . 


- I\| 2 A 'nk, k 


m 1,•••,p/2 


Equations (29*) and (29b) are a system of p 11 ne a r homogeneous 
equations for the determination of the A* and A* * For a non¬ 
trivial solution, the determinant of the coefficients must vanish, 
namely, 


1 1 

H 

1 

• 

• 

• 

• 

-l 

a _ -iw a -iw 

nl 1 nl 

a nt * iM 1 

a nl +U *l 

1 

-1 

see 

a , -i<*» 

nl 2 

• 

• 

• 

1 

0 . + iw, 
ni 1 

• • 

• 

°nl ***p /2 

(n + iw )i 
e no o 

(a _+ iw )t 

6 nl o 

(-One + H>) T 

V 

a . + iw 

nl 1 

a . + iw 

nl l 

-a . + iw 

ni 1 


# # e 


(30) 


- 0 


From this equation ve determine the 
*WviM requires that ve know the functions 


A , n ■ 1,2 ,••• 
n 


°nk - W 

These functions f (X ) are determined by solving for the roots 
of (25) in terms of ft . We can do this analytically in terms 
of algebraic functions for equations of fourth or lover degree. 
t’vh a would that problems could in principle be h a ndle d for 

D((iw) 2 ) of eighth or lover degrees. 


(31) 
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Equations (29c) and (29d) vill be found to be Identically 
satisfied. All constants A+. , k » 1, .. p/2-1, A^y, k - 1, 
p/2 can be solved for in terms of A+p /2 using equations ( 29 a) and 
( 29 b). -A^p /2 Is determined from the normalizing condition 



The case where the power spectral density Is a rational 
function of frequency Is a very Important one and covers a 
vide range of practical cases. 

b) Orthogonal Expansion of Processes with Markoff 
Autocorrelation Functions 


Suppose that the autocorrelation function of a stationary 
random acoustic noise process p(t) of mean zero is given by 


B (p(t) p(t +t)} - R(t) ■ A e 


-JJ|t 


(33) 


A ■ E(p (t)} ■ R(o), p » constant 


This function Is shown In Figure 3» 



Figure 3 . 


Graph of the Autocorrelation Function Ae 
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The spectral density corresponding to (33) 


S(w) 



e -p|t| - i*T 



2J$ 

2 2 

-(i »)% 


This function Is plotted in Figure 




Figure Spectral Density Function of a Process with Markoff 

Autocorrelation Function 


From (34]t, we see that the process under consideration has a 
rational spectral density function with 


*Ui" f) - a?* 

DUi») 2 ) - -(i *) 2 + @ 2 

The differential equation (19) becomes in this case 


2PA* n (t) 


|a| 2 (- 


♦ p 2 *J 

at 2 


Let 


« 2 . -a 2 + , 

“ l> ' 2 

n 


Then ( 36 ) can be written as 



(36) 

(37) 
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(38) 



» 2 t 

n n 


0 


“ p2< tV' 

' n' 


then (^ 6 ) has solution 


t„(t) 


iw t -iw t 

ae n + b e n , w > 0. 
n n n 



a n 5114 b n *** “^ttrary constants. Substituting ( 33 ) into 
the Integral equation ( 11 ), we find 


t I 

A j' e^t^t'Jdt' + aJ sP T t n (t , )4t' 


\\\ * (t), T - t-t 
n n 


Using ( 39 ) and carrying out the integrations glres 


e i« a t -pt -it^t -pt 

Aa - ■ * + Ab -■ + 

n n 


P + iw 


P - iw 


+ Aa 


{ p(t-T) + IWqT 
-p + iw 


+ Ab 


,P (t-T) - lUjjT - -iw_t 
^_ m 6 _ u 

-P - lM 


|x I 2 a e^ + |x I 2 b e'^ 

n n q' 

Ibis equation is identically satisfied in t if 

Aa Aa 

p _ n 1 , .2 

P + iw *p + iw n a n 

n n 

Ab Ab 

- 2 —_ 2 . K | 2 . 

p - iw -p - iw 'V n 

n n 
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-Aa 


Ab 


P + iw 


P - iw 


. (-0 + iw )T 

Aa e n 

n _ 

-P + iw 


Ab .-O + *“«>* 
n _ 

p + iw 


- 0 


The first two conditions in (40) are identically satisfied. Die 
second two are two linear homogeneous equations for a b . 

For a nontrivial solution, we require n n 


1 



P + iw 
n 



(-P + iw )T 
e n 

-p + iw 
n 


^-(P ♦ l« n )T 

P + iw 
n 



This reduces to 


1 1 

w cos - w T + p sin — w T » 0 
n 2 n 2 n 


u 

n 






Equation (1(2) determines the eigenvalues A . It Is clearly 
a complicated equation and must be solved numerically. 

From (40) 


b 

n 



P + Iw 
n 


a 

n 


(*3) 


Therefore 
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♦_(t) 


P+1U 


) 



<<*> 


P+1. ( “n c “ v + P « t) 


From the normalizing condition we find 


|a | 
n 


ft + in 


if |« n CO. « n t + P sin W Q tI 2 dt} 


-1 


(* 5 ) 


•the orthogonal expansion of p(t) now take, the fora 


P(t) 


I 


i\i« n 

rnz co * V + P v> X n 


(46) 


x n can only he determined when a realized value of p(t) 
is given. * 011 . In principle complete, the problem. 



Another method of solving the Integral equation (ll) much 
used by physicist. Is to expand the kernel R(t-V) Into a suit- 
able set of orthornormal functions. By this weene the problem 
can be reduced to the problem of finding the non-zero roots of 
an Infinite determinant. By Judicious choice of the expansion 
this process can be made tractable. (Ref. 5 .) 


d) Other Solutions 

1* Stationary Band Limited Qausslan White Noise of 
Mean Zero (Ref. 4 .) 


The expansion Interval in -• < t < « . The autocorrelation 
function for the process Is 


«/_\ S 2 N »in 2 ttBt 

r(t) • Cp > ~^r~ - 

and Is shown In Figure 5 , 
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< p > I 




Figure 5. Autocorrelation Function of Band Limited White Noise 

B ** w - is the bandwidth of the noise. spectral density 

appears as in Figure 6. 



Figure 6. Spectral Density of Band Limited White Noise 


Note that <Cp S B. The proper functions for the expansion of 
p(t) on the interval (•>) are found to be 


♦ (t) -«/3 

n 2irB(t - n/2B) 


m 
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{* ) is an orthonormal set but t Are not solutions of Eq. 11. 
n n 


The expects ion takes the form 


P(t) - 


T 


sin foB(t . n/2B) x 
2uB(t - n/2B) “ 


(*9) 


where 


-1 


X - 


atRR)' frit) 

iL 2rrB(t-n/2B) 


dt 


(50) 


2. Band Tilted Stationary formal White Boise on a 
finite Interval. (Ref.6.) 


Ihe correlation function is still of the form (Vf), but the 
expansion Is desired on the Interval (o,T). The proper functions 
and eigenvalues are 


♦ (t) 

n 


[*BT, 2/T(t-l)) 
on 



o < t < * 


X 

n 


H 

on 


fp 2 > KB 1)J 2 

on 



2J + 1 



R (c,y) and 3 (c, cos 0 ) are refpectively, the prolate 

ni nl 

spheroidal Bessel function and prolate spheroidal Legendre 
function as given by Morse and Fesbbech (Ref. 5). These 
functions have been tabulated. The expansion reads 
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X are Independent noxml variates of mean zero, 
n 


3. White Molae 

Correlation and Spectral density functions are 


R(t) - S q B(t) 
S(«) - S q 



S q Is a constant. The Integral equation In this case becomes 



T 

S B(t-t') t (V)dt' 


l\l* ♦,<*> 


or 


8 ♦ (t) - |X | t (t) 

on n n 


This can be satisfied by any orthonormal set of functions If 

we take X ■ vS for all n. The expansion becomes 

no 


p(t) -n/s ) t (t) X 

oZ_j n n 

T 

\ ■ JT / P(t) ♦n (t)dt 



The results for white noise are very Interesting and 
potentially wary useful. We can make the general statement 

that very broad band random signals can always be approximately 
expanded Into orthogonal series using any suitable set of 
orthonormal functions. 

17. Methods of Approximating Acoustic Boise Fields 
a) Optimal Approximation Techniques 

Once X and t (t) are determined from solving the Integral 
equation, $ can §e obtained when a typical recording p(t) is 
available. & We get a set of realized values for the X . The 

series ) X ♦ (t) ■ p(t) Is then an analytic representation of 

the realized function p(t). The sirens could be used to try 
to simulate this series or the typical recording. 
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In general, when vt hare an analytic expression for p(t), we 
need a synthesis procedure to approximate p(t) with the sirens. 
We might try to approximate p(t) Itself in the time domain by 
requiring the mean square difference between the sirens' output 

> p t (t> and p(t) to be a minimum. Oils requirement is achieved 
by minimizing 





with respect to the adjustable parameters of the sirens such as 
center frequencies, amplitude of outputs, percent modulation 
Let X be a parameter of the kth siren, then this synthesis 

procedure requires 



N is the number of sirens and L Is the number of adjustable 
parameters of the k^ siren. (j?6) leads to a system of equations 
for the optimizing parameters. Unfortunately, the system is 
nonlinear in the parameters in almost mil cases and it is difficult 
to solve even with digital computers. Here we try to optimally 
synthesize one realized value of p(t). Since p(t) Is a sample 
function from a random process, it may be better to approximate 
some statistical characteristics of the process. 

Other synthesis criteria might be: 

1. Minimize mean square difference between power spectral 
density of p(t) end that of y p k (t>. 

2. Optimize probability density of p(t). 

T N 

3. Minimize | E{f [p(t) -^P k (t)] 2 dt} 

® k-1 

with respect to the parameters of the sirens. We assume here that 
some of the parameters of sirens are random in character. 

b) Spectral Broadening of a Single-Tone Siren . 
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Spectral broadening of a single-tone siren can be achieved 
by various forms of modulation or by switching techniques. Let 
us consider first spectral broadening by modulation. 

1. Amplitude Modulation (Deterministic) 

A single-tone siren output may be represented by 


p(t) - A cos (« c t + f) (57) 

where A Is the amplitude of the output, u is the center fre¬ 
quency of the siren, and f is a phase angle. Such a siren output 
as (57) contains only one spectral component. Suppose by some 
means we are able to vary A as a function of time. In particular, 
let 


A(t) • A (1 ♦ A cos w t) 
o m 


(58) 


A , A, and u 

o m 


are constants 


Then (57) reads 

p(t) > 


tA 


4 


tv •' 

ur. 


Of) Ouj Met Y - 


t 

tJc t jv-v 


A (1 + A cos w t) cos (u t + f) 
o me 


(59) 


or 

p(t) • A cos(w t+ f) + — A A cos [(« + w )t + p) + 
o c 2 o cm 

x (60) 

+ - A A cos [(w - id i ) t + f] 

20 cm 

By this modulation, we see that p(t) now has three spectral 
components. It is easy to see that If 


1 

P<t) " A oI (1 ♦ \ c °* COa («c t + *) ( 61 .) 

k-1 

p(t) will have 21 + 1 spectral components. Equation (6l) includes 
a wide variety of amplitude modulation such as arbitrary periodic 
modulation and almost periodic modulation. He note, however, that 
as long as A(t) is a deterministic function of time the spectrum 
Is discrete. 


22 






2. Frequency Modulation (Deterministic) 


The instantaneous a n g ul a r frequency is defined as 

“<*> • JJ **» (fie) 

vhere the phase angle f Is now considered a function of tine. 

If me vary the frequency by a single tone, -On siren output 
takes the font 


p(t) ■ A cos[w t + sin w tl 

o CM m 


( 63 ) 


A ind m ax# constants 


It is well known that t hi s can be written in the fore 


p(t) ■ ~ A COS M t + i 


| f 1 + ( -1 ) 1 J n (“ ) (00s (M c t ♦ nw^)t 


+ coa (m - DM )t) 

e a 


«*) 


from which we see that there are an Infinite uuaber of discrete 
spectral components. Is the Bessel function of the first 

kind of order n. general appearance of the Intensity spectra 
of p(t) for small. Intermediate, and large values of A /m are 
shown in Figure 7 with their envelopes. * 




















For arbitrary periodic or almost periodic frequency modulation, 
the siren output is given by 



I 

A cos {« t ♦ ) 

o c Z_i 

k-1 



sin w 


mk 




A and w are constants, 
dt mk 

The number of sidebands increases tremendously as H increases, 
but the intensity apectrua is still discrete. Only for an 
aperiodic modulating save can we obtain a continuous intensity 
spectrum or spectral density. 


3.) Amplitude Modulation (Random) 

Continuous spectral broadening of a single-tone siren can be 
achieved through random modulation. Consider 

p(t) « A [1 + A P w (t)] cos w t (66) 

O n c 

where P (t) is a stationary normal random process and A , A, and 

w are constants. Bn spectral density of p(t) for a modulating 
noise with spectral density function given by 

s (») - s «*> 

M O 

S u m constants 

o b 


is shown in Figure 8 for various modulation Indices. 

It is seen that spectral boradening can be quite considerable. 
Rote that finite energy is in the carrier for all modulation indices. 

shown. 


4.) Frequency Modulation (Random) 
In this case we have 


p(t) - A cos [« t + P (t)] 

O C A 

where P (t) is a stationary normal random process. A q and « c are 

constant. The spectrum of p(t) for P^t) taken to have spectral 
density 


( 68 ) 











3 m M 


k <*> \ 

i + « 2 /b 2 

' V 


V < > 6 « 

tin ( Q A + *>Z ) (^ 9 ) 

£«/ 


i* shown In Figure 9 for various conditions. B is the band- 
vidth of the noise. <1^>is mean square value *of modulating 
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Figure 8, Spoctrtl Density of Noise Modulated 

Single-Tone Siren 
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If n < < l, the spectral density is very narrow shout *» c 
(little broadening). For very large values of p, the 
spectral density i® bread about w . 

5.) Spectral Broadening by Random Switching 

Consider a stationary noise pressure of the type shown in 
Figure 10. 



The zero crossings occur at purely ra ndo m instants t^,t^,.... 
Such an output would be produced (except for a carrier ten* 
cos w t) by an ideal siren that vas being switched on and off 
at purely random times. The spectral density function of the 
pressure in Figure 10 is 


T 

S(w) ■ ~ 
H ir 


1 ♦ 


€*f) 


(70) 


vhen t^-t* is a Poisson distributed random variable. 

If ve use l(t) as an amplitude modulation function of a single* 

tone siren output; te have the sound pressure 


P(t) - H(t) cos [mt ♦ f] 


jj ^ random phase uniformly distributed between o and 
autocorrelation function of p(t) is 


(71) 

























R - E{l(t) cos l* c t + f] f B (t + t) cos [« (t + t) 4 f]) (72) 

c 

1 r 2 * 

R p • E{H(t) H(t + t)} — J cos [« t + ?] cos [m (t + t) +9 ] df 

o c 

V t) - 2 "» (t) co * V 

The spectral density is obtained by use of the w *n«. r . K-j n ^v .j 
theorem. We find 

s p ( w ) m J R p (*) e^ihr 

® ( w ) i r f )xi -(« )xl 

P * c + e c' ) RjCTjdr (73) 

s p ( w ) ■ 4 Cs H (»-« c ) 4 S^m + « c )} 
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V. Comclusloap and Beconmendatlona for Farther Study 


The energy of a single-tone siren can be spread over a fairly 
large frequency band about the center frequency of the siren by 
various sorts of modulation. The most effective say of spreading 
the sound energy over the frequency band Is by random amplitude 
or frequency modulation • 


A given acoustic noise field may be approximated (in spectral 
density) by centering the sirens' center frequency at those which 
constitute the given signal or by dividing the spectral density 
of the given signal Into equal energy bands and assigning an energy 
band to each siren. The approximation can be improved by amplitude 
«nH frequency modulating each siren. Approximation of a given 
acoustic noise field in some optimal way le a d s to the solution 
of a system of non-linear equations. 


Further work may be carried out to find other solutions to 
the integral equation occurring in the Karhunen-Iosrve expansion 
of the random noise field. It would be particularly desirable 
to find some exact solution so that the proper expansions could 
be compared with a Fourier expansion of the field. Another 
direction in which a continued investigation might bear fruit is 
that of trying to devise methods of solving the nonlinear equations 
arising in trying to optimally approximate bosk characteristic of 
the given acoustic noise field with the sirens. Still another 
area worth looking into further is spectral br o ad enin g by ra ndo m ! 
zation of the sirens' parameters. 
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